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Using bivariate and multivariate morphometric analyses, effects of paedomorphosis 
on intrapopulation ontogenetic changes in size and shape of three European newt 
species (Triturus vulgaris, T. alpestris and T. carnifex) are studied. Samples of larvae, efts 
and adults of two conspecific populations of each species were considered, one with 
generalized life cycle, the other with paedomorphosis (in a T. carnifex population 
larvae had prolonged growth and metamorphosed during the following year). It was 
found that existence of an alternative life cycle rendered more complex ontogenetic 
picture of allometric relations of individual morphometric features, particularily in 
T. vulgaris. In this species paedomorphic individuals followed ontogenetic trajectories 
of larvae, while metamorphic adults deviated considerably from this trajectories. This 
was much less significant in T. alpestris, another species with high incidence of paedo- 
morphosis, but also in T. carnifex, a species with less plastic life-cycle than both 
previous species. As in the American newts (Harris 1989), differences between non- 
transforming and transforming individuals in the genus Triturus in ontogenetic paths 
of size and shape changes presumably developed relatively late in ontogeny. 

Ana Djorovic, Milos L. Kalezic, Institute of Zoology, Faculty of Biology, Studentski 
trg 16, 11000 Beograd, Yugoslavia. 


Introduction 


It is the commonly held view that the heterochrony, changes in the timing of developmental events 
during ontogeny, is among the most prominent mechanisms for phenotypic changes (e.g. McKinney 
1988). Allometry, the study of size/shape relations, provides data for assessing heterohronic morpho- 
logical events (Alberch et al. 1979). If paedomorphosis, sexual maturation of individuals retaining 
morphology that resembles subadult stages of closely related populations or species (Gould 1977), is 
to be studied, ontogenetic allometry is relevant. Ontogenetic allometry simply refers to changes in size 
and shape relationships in individuals with time (Shea 1984, Lande 1985). 


Paedomorphosis is frequent in amphibians, especially in Caudata where it occurs due to change in 
timing of metamorphosis. Metamorphosis consistently occurs in generalized and presumably primitive 
form of urodele life-history; from eggs laid into aquatic environment larvae are developed and 
metamorphose into immature individuals (juveniles or efts) who remain mainly terrestrial before 
reaching sexual maturity. In many urodele species, however, there are populations with individuals 
following another ontogenetic pathway. Larvae here have prolonged growth, attain sexual maturity 
and reproduce retaining larval morphology. Immature eft stage is excluded in this path, while adults 
can have abandoned or delayed metamorphosis. These urodele species exhibit two adult phenotypes 


in natural populations. Paedomorphs are adults having such larval characteristics as: external gills and 
associated hyobranchial structures, an enlarged tail fins, skin often with Leydig cells and larval texture, 
and characteristics of larval skull. Metamorphs lack all these traits in morphology. 

There is considerable interspecific variation among urodeles with paedomorphosis in the degree of 
retention of larval characteristics. Reilly (1987) found, studying metamorphosis of skull and hyo- 
branchial apparatus, that paedomorphosis in Ambystoma talpoideum has profound effects on these 
osteological structures, while paedomorphosis in salamandrids (Salamandridae; Notophthalmus viri- 
descens, Taricha granulosa and Triturus vulgaris) is mainly limited to the retention of gill structures to 
varying degrees. Also, it seems that urodeles species with paedomorphosis differ in the stage of somatic 
truncation along an ancestral size-shape trajectory (Shaffer 1984, Harris 1989). 

Paedomorphosis is rather common in many species of the genus Triturus. The smooth newt (Triturus 
vulgaris) and the alpine newt (Triturus alpestris) show a high degree of paedomorphosis incidence, 
particularly in the submediterranean area of the Balkans (Dzukic et al. 1990). In the crested newt 
(Triturus cristatus superspecies) almost an exclusive life-history pathway leads to obligate metamorpho- 
sis; i.e. paedomorphosis is extremly rare. So far, there has been no study on paedomorphosis related 
ontogenetic changes in size and shape of European newts (genus Triturus). The main goal of this paper 
is to compare such intrapopulation changes in three species, T. vulgaris, T. alpestris, and T. carnifex. 
Populations with paedomorphosis were studied, as well as conspecific populations with obligate 
metamorphosis considered here as control groups. The null hypothesis in this study was that the 
allometric relations of paedomorphs should resemble more those of larval and juvenile individuals 
than those of metamorphs ones. 


Materials and methods 


Population Samples 


Paedomorphic population of Triturus vulgaris from the pond of Velika Osjecenica (Vilusi, Montene- 
gro), further designated as V1, was chosen because of high percentage of paedomorphic individuals 
(Tucic at al. 1985). Morphometric variables were scored on 205 larval and post-larval specimens 
collected from this pond. The smooth newt population from pond Cesta (Vrsacki breg, Serbia) was a 
control population (population referred to as V2), as no evidence of paedomorphic individuals was 
found in the years of its inspection. When larvae from this pond were sampled (July), some were 
sacrifised and preserved as larval sample, while others were grown under laboratory conditions in 
order to provide eft sample. In the end of August the majority of larvae metamorphosed; efts were 
grown for two weeks in terrarium and then preserved. The total population sample contained 111 in- 
dividuls. 

The alpine newt population from Bukumirsko Lake (Montenegro; population referred to as A1) was 
selected as paedomorphic population. In this lake non-transformed individuals predominated (Ra- 
dovanovic 1951). Neither larvae nor efts were found in this lake. Efts were sampled in a nearby 
temporary lake (Crno Lake) separated from Bukumirsko Lake only by a low moraine which presents 
no barrier to intensive gene flow between the newts of these lakes. We failed to collect larvae from this 
lake, too. The alpine newt paedomorphic population sample consisted of 169 individuals. Samples of 
a control alpine newt population (referred to as A2) were taken from the locality of Sirokar (Montene- 
gro, Mt. Prokletije). From here we took a total of 166 newt specimens. 

The individuals from Lokanj pond (Montenegro), referred to as C1, have plastic life-history, includ- 
ing the strategy of paedomorphosis (Kalezic et al. 1994). We have sampled here, besides efts and adults, 
overwintering larvae as well. The total population sample estimated 118 newts. The individuals of 
crested newts from nearby Ceklin pond are characterized by obligate metamorphosis. This population 
is a control one in our study and is referred to as C2. From this population morphometric characteristics 
of 179 individuals were measured. 


Measurements 


To demonstrate the size and shape changes during newts ontogenesis, the following set of 10 mor- 
phometric variables were measured to the nearest 0.1 mm using dial calipers or plastic ruler as 
appropriate: L: total length, Lep: body length (measured from the snout to the front edge of cloaca), Led: 
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tail length, Lsv: snout-vent length (measured from the snout to the posterior edge of cloaca), Ltc: head 
width, Lc: head length (measured from the snout to the angle of the jaw), Pa: forelimb length, 
Pp: hindlimb length, D: distance between fore- and hindlimbs, and Lh: maximum tail height (measured 
at the base of the tail from dorsal to ventral edges). 

Before they were mesured, the individuals used in this study had been preserved in ethanol for 
different period of time (from a couple of months to 7 years). In scoring morphometric traits values and 
their on-going statistical treatment, we may have dealt with a source of potential errors. We are aware 
that occasional imprecision in measuring might occur due to some morphological distortion which can 
affect newts during their preservation (Lee 1982, Verrell 1985). The extent of this measurer bias in our 
study was not determined. But, as individuals of each ontogenetic stage spent approximately the same 
time in conservation solution, and as it might be supposed that tissue shrinkage due to preservation 
was not age and paedomorphs/metamorphs dependent, we may assume that this potential error did 
not introduce any systematic bias into analysis. 


Statistical analysis 


Since we have treated ontogenetic stages (larvae, efts, adults) of the same population as different 
categories, obtained sets of measurements could be referred here to as “cross-sectional data” (Cock 
1966). These data yield information on the average growth of individual groups. Due to marked sexual 
size differences in newts (Kalezic et al. 1992), sexes of paedomorphs and metamorphs were separately 
statistically treated. 

Univariate traits were analyzed using one-way analysis of variance (ANOVA) when comparing the 
significance in difference between trait variability, and two-way analysis of variance when evaluating 
the significance of paedomorphosis and sex effects on morphometric character diversity. The used 
measure of bivariate allometry is the regression coefficient of logarithmically transformed data. In this 
analysis, Lsv was used as an independent size variable. The differences among bivariate allometric 
coefficients of morphometric characteristics for analyzed forms were tested by Tukey q-test for multiple 
slopes comparison (Zar 1984). 

Patterns of morphological variation were examined by the multivariate methods of principal com- 
ponent analysis (PCA) done on pooled all measurments data. The first component (PC1) of this 
analysis, which is highly positively correlated with original data, was used as a latent size variable, 
while the second principal component (PC2) measures organism shape independent of size (Shea 1985, 
Bookstein et al. 1985 and references therein). The analysis was performed on the variance-covariance 
matrix (more appropriate when one deals with morphometric data than correlation matrix) of log- 
transformed variables (done in order to meet the assumption of homoscedascity; Zar 1974). 

A test of isometry was performed by comparing PCA eigenvectors to those of an isometric vector 
(chi-square test; see Jolicoeur 1963). To measure relationship between bivariate and multivariate 
allometry, Spearman’s rank correlations were calculated between the eigenvector loadings for the first 
component and allometric loadings of each morphometric characters. 

Ontogenetic trajectories were obtained for each ontogenetic stage (larvae, juveniles and adults) by 
linear regression analysis of PC1 scores on PC2 scores. A rate of change in shape relative to size was 
estimated by the slope of the regression lines. 


Results 


Morphometric traits variability 


The mean values (+one standard error) of 10 morphometric characters for larvae, efts, and adults 
(females and males were given separately within paedomorphs and metamorphs groups) of examined 
populations were given in Appendix for T. vulgaris (Tab. 1), T. alpestris (Tab. 2), and T. carnifex (Tab. 3). 
In the smooth and crested newt population samples efts were significantly larger than larvae individ- 
uals for most characters (ANOVA test), while such differences in the alpine newt were much less 
pronounced. Among the adults, in the smooth newt paedomorphic individuals were significantly 
smaller than metamorphic ones for both sex; such differences within alpine newt population were 
much less apparent. 


The influence of paedomorphosis effect and sex effect on morphometric variability in the smooth 
newt appeared to be statistically highly significant for all characters, except for tail height (two-way 
ANOVA test). In the alpine newt such influences were more selective. Intersex differences in examined | 
morphometric characters were much more pronounced in the crested newt population with obligate 
metamorphosis than in the population with more plastic life-history traits (ANOVA test). 


Bivariate allometry 


In the smooth newt population with paedomorphosis (V1) larvae retreated significantly from isom- 
etry for many characters, following with metamorphs and efts in decreasing order (Tab. 1). In paedo- 
morphic individuals statistically significant coefficients of bivariate allometry were only in some 
characters (head dimensions and tail height). In the control smooth newt population (V2) all measured 
morphometric characters of males had isometrical growth, following by larvae where only tail height 
had allometric growth. Among efts and females allometry appeared in more characteristics (mainly | 
head dimensions). 

In the alpine newt population with paedomorphosis (A1 ), a number of characters of non-transformed 
females and some characters of transformed females, expressed statistically significant allometric 
growth (Tab. 2). Morphometric characters in males (paedomorphic and metamorphic) were almost 
excusively characterized by isometry. However, in the population of this species with obligate meta- 
morphosis (A2), most characters of larvae and males were with a significant bivariate allometry. In both 
crested newt populations, the number of significant character allometry decreased from larvae, 
throught efts to females and males (Tab. 3). 


Tab. 1. The significancy of allometric coefficients - testing the slope against a standard 1 - for all analysed form in 
smooth newt populations. n.s.: isometry; statistically significant allometry *P<0.05, **P<0.001, ***P<0.001. See text for 
character and population designation. 


Vi V2 
metamorphic paedomorphic metamorphic 

larvae efts ee 3d el 3d larvae efts IQ) 3d 
IL; at ns. n.s. n.s. n.s. n.s. n.s. n.s. n.s. n.s. 
Lep bs = ides ae n.s. n.s. n.s. n.s. ie n.s. 
Led = n.s. n.s. n.s. n.s. m.s: n.s. n.s. n.s. n.s. 
Lte So ee $ n.s. EE n.s. n.s. DoE aks n.s. 
Le bt * A * ae ++ n.s. ns: kan mess 
Pa bas n.s. b isa n.s. n.s. n.s. n.s. n.s. n.s. 
Pp EEE ii a be n.s. n.s. n.s. r n.s. n.s. 
D n.s n.s n.s n.s n.s n.s n.s. n.s. n.s. Dsi 
Lh biy $ bi Soh Er n.s. 3 n.s. n.s. n.s. 


Tab. 2. The significancy of allometric coefficients - testing the slope against a standard 1 - for all analysed forms in 
alpine newt populations. n.s.: isometry; statistically significant allometry *P<0.05, *P<0.001, ***P<0.001. See text for 
character and population designation. 


Al A2 
metamorphic paedomorphic metamorphic 

efts 22 3d ee 3d ' larvae efts ee feo) 
E n.s. n.s. = aiii n.s. g n.s. n.s. x 
Lep n.s. n.s. n.s. = n.s. ns. n.s. ns. = 
Led n.s. n.s. n.s. = n.s. = n.s. n.s. F 
Ete axe n.s. n.s. n.s. n.s. sed n.s. n.s. n.s. 
IL n.s. xe n.s. n.s. n.s. os n.s Š Si 
Pa ns + T.S * n.s ae + He + 
Pp $ n.s. n.s. iid nS. z n.s. Š n.s. 
D n.s. n.s. n.s. . n.s. n.s. n.s. n.s. n.s. n.s. 
Lh n.s. ` y * ns. n.s. $ Aa n.s. n.s. 
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Between-ontogenetic-stages pairwise comparisons in allometric coefficients were presented in Tab. 4 
in the form of the number of statistically significant differences of nine morphometric variables versus 
snout-vent length for each population. In both populations of the smooth newt and control population 
of the alpine newt (A2) there were no statistically significant difference between bivariate coefficients 
of larvae and efts, but it was not the case in the crested newt populations. Among larvae and adults, 
most differences existed between larvae and metamorphosed males, especially in the smooth newt 
population with paedomorphosis (V1). In the population with paedomorphic individuals, efts and 
paedomorphic females differed in two of nine comparisons. In these populations there was no differ- 
ence in regression coefficients between efts and metamorphosed females. The difterence between efts 
and metamorphosed females existed in both populations of the crested newt. Efts and males differed 
only in the populations V1 and C1 (Tab. 4). Statistically significant differences of bivariate allometric 
coefficients between paedomorphs and metamorphs of the same sex were found in the smooth newt 
and in the alpine newt only for one character comparison (Tab. 4). Intersex differences within meta- 
morphs existed in the smooth newt population V1, and in populations of the crested newt, while 
paedomorphic females and males differed in the population A1 for two bivariate relationships. Taking 


Tab. 3. The significancy of allometric coefficients - testing the slope against a standard 1 - for all analysed forms in 
crested newt populations. n.s.: isometry; statistically significant allometry *P<0.05, **P<0.001, ***P<0.001. See text for 
character and population designation. 


ell (G2 

larvae efts og fee) larvae efts 0o 3d 
L x * * 3 n.s. n.s n.s n.s 
Lep ns. n.s. n.s. ns. n.s. n.s. n.s. n.s. 
Led k pE b n.s. ns. n.s. n.s. $ 
Ltc ns. nS: n.s. n.s. et te ns. n.s. 
Le big n.s ms: n.s aht Ee n.s. n.s. 
Pa n.s bi n.s. n.s i n.s sipi n.s 
Pp e n.s. sabes n.s bz DS: te n.s 
D x TES: n.s n.s n.s x n.s. n.s. 
Lh n.s. n.s. ns. * n.s. a n.s. n.s. 


Tab. 4. Number of statistically significant differences between bivariate allometric coefficient (q-test) of nine 
morphometric characters and Lsv as independent variable of larvae, efts, paedomorphs and metamorphs (sexes were 
separated) in six studied populations. Abbreviations: L: larvae, E: efts, FM: metamorphic females, FP: paedomorphic 
females, MM: metamorphic males, MP: paedomorphic males. df: degree of freedom. See text for population designation. 


V1 V2 A1 A2 C1 (A 
df=297 df=103 df=168 df=158 df=110 df=171 
L/E 0 0 if 0 il 2 
L/FM 1 0 / 0 0 0 
L/FP 0 / / / / / 
L/MM 4 1 if 0 0 it 
L/MP 1 it ii / if / 
E/FM 0 0 0 0 2 1 
E/FP 2 / 2 / / / 
E/MM 2 0 0 0 ‘l 0 
E/MP 2 / 0 if / / 
FM/FP 1 if 0 if / i, 
FM/MM 1 0 0 0 1 A 
FM/MP 1 i il / / if 
FP/MM 1 / 1 j / / 
FP/MP 0 / 2 / / / 
MM/MP 0 / 0 / / if 
Total 16 1 6 0 5 6 
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Tab. 5. Character eigenvectors and eigenvalues derived from pooled data of larvae, efts, paedomorphs and metamorphs 
of the six populations. The variables used in the analysis were the logarithms of 10 morphometric characteristics. See 
text for character designation. 


Lh 


Eigenvalues 


% of variance 


BEI 


0.31 
0.31 
0.31 
0.31 
0.28 
0.23 
033 
0.34 
0.36 
0.35 


0.3825 
(93.42) 


Eigenvectors 
REZ PC3 
-0.06 0.31 
-0.16 -0.03 
0.03 0.66 
-0.15 -0.02 
-0.10 -0.58 
-0.16 -0.14 
-0.13 -0.22 
-0.14 -0.12 
-0.15 0.18 
0.92 -0.13 
0.0195 0.0031 
(4.77) (0.76) 


Tab. 6. Chi-square values obtained from test of isometry (comparing PCA eigenvectors based on covariance matrix 
with an isometric vector -arc cos 1/3; total df=9). The variables used in the analysis were the logarithms of 
10 morphometric characteristics. See text for character designation. P: significance of differences from isometry. 


*P<0.05; ***P<0.001 


Larvae 
Efts 

22 (M) 
3d (M) 
22 (P) 
3d (P) 


V1 


2.94 
21.06* 
15.14 

9.69 

3.09 

4.97 


V2 


54.88*** 
392 
4.52 
4.74 

/ 
/ 


Al 


/ 
17.43* 
54.58*** 
ISS 
58.64*** 
12.67 


A2 


1.66 
2.34 
5950 
11.51 
if 
/ 


C1 C2 
5:99 4.88 
8.39 312 
1.47 2.84 
575 3.03 
/ / 
/ if 


Tab. 7. Spearman’s rank correlation between results of bivariate and multivariate allometry. M: metamorphic 
individuals. P: paedomorphic individuals. 


Larvae 
Efts 

99 (M) 
36d (M) 
22 (P) 
3d (P) 


Vil 


0.720* 
0.500 
0962 
0.950*** 
0.946** 
0.427 


V2 


0.418 
0.284 
0.83377 
0.828** 
/ 
/ 


Al 
/ 
0979 
0.883** 
0.770* 
0.200 
0.217 


A2 


0.300 
0.736* 


-0.017 


0.703* 
7 
/ 


Cil C2 
0.633* 0.833** 
0.845** 071 
0.600* 09662 
0.946*** 0.9675 

/ / 
/ / 


Tab. 8. Means and standard errors of PC1 and PC2 scores for analysed forms in the smooth newt populations. See 
text for population designation. 


Larvae 
Efts 

22 (M) 
3d (M) 
TSE) 
3d (P) 
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Rel 


3.06 + 0.06 
2.92 + 0.05 
4.13 + 0.07 
4.16 + 0.09 
4.04 + 0.06 
3.10 + 0.07 


V1 


PC2 


-0.31 + 0.03 
-0.71 + 0.02 
-0.46 + 0.05 
-0.39 + 0.04 
-0.37. + 0.04 
-0.38 + 0.04 


IRCI 


2.67 + 0.10 
2.62 + 0.08 
3.89 + 0.08 
4.05 + 0.11 
/ 
/ 


V2 
PEZ 


-0.40 + 0.10 
-0.70 + 0.03 
-0.57 + 0.05 
-0.36 + 0.05 
/ 
/ 


into account the total number of statistically significant differences between bivariate allometric coeffi- 
cients of all nine morphometric characters in larvae, efts and adults comparisons, some trends become 
apparent (Tab. 4). Populations with paedomorphosis (T. vulgaris and T. alpestris) were much more 
different from conspecific populations with obligate metamorphosis than conspecific populations with 
much less differences in life-history pathways (T. carnifex). 


Multivariate allometry and ontogenetic trajectories 


All morphometric variables were positively correlated with first principal component (PC1) which 
explained more than 93 % of the total variability in six newt populations (Tab. 5). With the principal 
component two (PC2) only tail height (Lh) and tail length (Lcd) were positively correlated. PC2 
component explained almost 5 % of total variability. Next principal component (PC3) took less than 1 % 
of total variation, and therefore was not further considered. 

The results of assessing multivariate allometry on separate PCA for each forms of analyzed popu- 
lations were given in tab. 6. The statistically significant differences from isometrical growth were found 
in paedomorphic population of the alpine newt for efts and both groups of adults, expect for paedo- 
morphic males. In the conspecific control population (A1) such a difference was found only for females. 
In paedomorphic population of the smooth newt a difference from isometry was found only for efts. 
All forms in populations of the crested newt had a multivariate isometrical growth. 

We also tested the correlation (Spearman’s rank correlation) between results of bivariate and mul- 
tivariate allometry (Tab. 7). These allometries were significantly correlated for metamorphic males and 
females (except in A2 population) in all analysed populations. Non-transformed newts, however, had 
decopuled bivariate and multivariate allometries, except for pedomorphic smooth newt females. When 
preadult stages were considered, only in smooth newt populations (efts) and in populations with 
obligate metamorphosis of the smooth newt and alpine newt (larvae), these allometries were not 
correlated. 

Means and standard errors of PC2 and PC1 scores for analyzed forms (larvae, efts and adults) were 
given for the smooth newt (Tab. 8), the alpine newt (Tab. 9) and for the crested newt (Tab. 10). 
Concomitant bivariate plots of PC1 and PC2 scores and regressions for these formes were presented in 
Figures 1 and 2 for the populations with paedomorphosis only. The most informative way to compare 
obtained ontogenetic trajectories of examined newt populations was to test significance of differences 
of the slopes of regression lines of PC1/PC2 scores for larvae, effts, paedomorphs and metamorphs 
within populations (Tab. 11). Some trends in these comparisons were obvious. Larvae did not differ 
from efts in the rate of change in shape relative to size in any studied newt species. Larval and 


Tab. 9. Means and standard errors of PC1 and PC2 scores for analysed forms in alpine newt populations. See text for 
population designation. 


Al A2 
Pcl PC2 PC1 PE 
Larvae / / 3.00 + 0.06 -0.31 + 0.03 
Efts 3.06 + 0.08 -0.71 + 0.03 2.85 + 0.06 -0.71 + 0.02 
22 (M) 4.22 + 0.08 -0.65 + 0.09 4,23 + 0.07 -0.65 + 0.08 
3d (M) 4.10 + 0.10 -0.61 + 0.07 4.10 + 0.06 -0.57 + 0.04 
Bene) 4,22 + 0.08 -0.56 + 0.09 / it 
go (P) 4.17 + 0.09 -0.51 + 0.07 i / 


Tab. 10. Means and standard errors of PC1 and PC2 scores for analysed forms in crested newt populations. See text 
for population designation. 


Ci @2 
PCI Pee RET EG? 
Larvae 4.22 + 0.05 -0.51 + 0.04 3.88 + 0.09 -0.49 + 0.04 
Efts 4.28 + 0.07 -0.60 + 0.03 4.02 + 0.10 -0.62 + 0.03 
22 (M) 4.76 + 0.10 -0.65 + 0.04 4.87 + 0.12 -0.61 + 0.05 
3d (M) 4.74 + 0.09 -0.64 + 0.04 4.74+0.11 -0.62 + 0.03 
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Larvae 
Efts 

P Females 
P Males 
M Females 


M Males 


o 1 2 3 4 5 6 
PC1 


Fig. 1. Bivariate plot of PC1 and PC2 scores and regressions for larvae, efts, paedomorphs and metamorphs (sexes 
were separated) for the smooth newt population (V1). 


paedomorphs’ (both sexes) ontogenetic trajectories were not significantly different in the smooth newt; 
larvae, efts, and paedomorphic females have almost parallel regression lines (Fig. 1). Ontogenetic 
trajectories of the metamorphs and larvae of the same population showed, however, significant 
differences. Efts in this population had quite different trajectories from metamorphs, but not from 
paedomorphs. Within-adult-sex comparisons revealed only significant differences between paedomor- 
phic and metamorphic females, but not between males. In the intersex comparisons paedomorphic 
females were apparently different from metamorphic males. 


Tab. 11. Significancy of differences between regression coefficients (q-test) of PC1/PC2 scores of larvae, efts, 
paedomorphs and metamorphs (sexes were separated) in six studied populations. Abbreviations: L: larvae, E: efts, 
FM: metamorphic females, FP: paedomorphic females, MM: metamorphic males, MP: paedomorphic males. n.s.: non- 
significant, *P <0.05, ** P<0.01, *** P<0.001. df: degree of freedom. See text for population designation. 


V1 V2 Al A2 C1 C2 
df=297 df=103 df=168 df=158 df=110 df=171 

L/E n.s n.s. / ns. n.s. n.s. 
L/FM bit n.s i. ae n.s. n.s. 
L/FP n.s i / / / if 
L/MM bas ns. / ns. n.s. n.s. 
L/MP n.s. / / / Ht / 
E/FM haii ns. n.s. aiza n.s. n.s. 
E/FP n.s. / n.s. / f it 
E/MM aks n.s. n.s. n.s. n.s. n.s. 
E/MP n.s. if n.s. 
FM/FP Ei / ns. it 
FM/MM n.s. n.s. n.s. n.s. n.s. n.s. 
FM/MP mes 7 n.s j / / 
FP/MM t if n.s ii if / 
FP/MP n.s. / nS if / / 
MM/MP , Te ii n.s / / / 
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Fig. 2. Bivariate plot of PC1 and PC2 scores and regressions for larvae, efts, paedomorphs and metamorphs (sexes 
were separated) for the alpine newt population (A1). 


Paedomorphic alpine newt population (A1), however, did not show any of such differences; there 
were no statistically significant differences in size and shape changes among analyzed stages, which 
was not the case with the control smooth newt population (V2) and with both crested newt populations 
(Tab. 11). In the alpine newt population with obligate metamorphosis (A2) the regression coefficient of 
females differed from that of larvae and efts (Tab. 11). 


Discussion 


A number of populations of the smooth newt (Triturus vulgaris) and the alpine newt (Triturus 
alpestris) from the submediterraneans of the Balkans have individuals with polymorphic phenotypes 
(see Dzukic et al. 1990 for review). Paedomorphs in these populations often outnumber syntopic 
metamorphs. Even a subspecific level has been attached to some of these alpine newt populations, 
including the alpine newt population with paedomorphosis studied here. This was based on phenotype 
differences of its paedomorphic individuals, including morphometric characters (head dimensions), in 
comparisons with the conspecifics of the nominotypical form which had obligate metamorphosis only. 

In the crested newt (Triturus cristatus superspecies) there are populations in which numerous larvae 
have prolonged growth and metamorphose into immature juveniles during the following year. How- 
ever, a few of the larvae with prolonged growth attain sexual maturity without metamorphosis, thus 
excluding a juvenile stage. But these individuals, contrary to other European newt species with 
paedomorphosis, do not stay long as paedomorphs. Up to now, in only one population of Triturus 
carnifex from Montenegro a few sexually mature newts with larval morphology have been found 
(Kalezic et al. 1994) - this population is included in this study (C1). 

We found in this study that larvae did not differ from efts in the rate of change in shape relative to 
size in any studied newt species. This is another confirmation that the main size and shape changes 
occurred after metamorphosis, till newts attain sexual maturity (Kalezic et al. in prep.). This was 
established to be about three years for the newts of all three species studied here from the Balkans 
(Kalezic et al. in prep., Cvetkovic et al. in prep.), and for the crested newt (T. cristatus) from France 
(Francillon-Vieillot et al. 1990). During this period newts mostly reach their adult size and shape which 
do not change significantly afterwards. But what happened when metamorphosis is omitted from life 
cycle? 


W 
N 
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In the smooth newt, paedomorphic individuals followed ontogenetic trajectories of larvae from the 
same population, while metamorphic ones deviated from this trajectories considerably. The morpho- 
logical null hypothesis postulated in this study that the size and shape of non-transforming newts was 
accounted for by an extention of the larval growth trajectries happened to be confirmed. The situation 
in the alpine newt population with paedomorphosis is much less clear probably because of inaccessi- 
bility of larvae from the population in our study. In this species trajectories of paedomorphic newts 
depicted no difference in comparison with trajectories of efts and, most likely, they would not do that 
of larvae too. But, metamorphs of the population did not deviate from this trajectories. Interestingly, 
the same happened with adults in other newt populations with obligate metamorphosis studied here. 
If any exception existed, as was the case in such a population of the alpine newt, females were those 
to deviate significantly from common population ontogenetic size-shape trajectories (Tab. 10). 
A general conclusion could be that paedomorphosis in some newt species (T. vulgaris in our study) has 
influence on ontogenetic trajectories making somehow metamorphic population counterparts different 
in this sense. In other cases, paedomorphosis (T. alpestris in our study) and more plastic life-history, i.e. 
prolonged larval period and omitted eft stage (T. carnifex), have no such effect. More plastic life-history 
had, however, an effect on interpopulation level; juveniles and adult males, the only sex with paedo- 
morphic individuals from Lokanj and Ceklin ponds, appeared to be statistically significantly different 
in the ontogeny of size and shape changes (Kalezic et al. 1994). 

Our results of bivariant analysis of shape changes during ontogeny are mainly in the line with the 
results of multivariate allometry. It was found that the existence of an alternative life cycle rendered 
more complex ontogenetic picture of allometric relations of individual morphometric features. The 
changes in allometric relations of paedomorphic and metamorphic individuals were particularly 
evident in the smooth newt. In this species paedomorphs more than metamorphs, followed the path 
of larval allometric coefficients. It seems also that paedomorphosis decoupled bivariate and multivar- 
iate allometries as the correspondence between these allometries at the adult stage was much more 
pronounced in transformed individuals than in syntopic paedomorphic individuals. 

It should be pointed out here that the sample of larvae taken from the population with paedomor- 
phosis in our study consists of individuals which would have metamorphosed, as well as of those 
which would have prolonged growth and reached sexual maturity without complete metamorphosis. 
There was no way of separating these two larvae groups. Due to this fact, as well as to using cross- 
sectional data instead of longitudinal ones, we were prevented from getting a clearer picture of onto- 
genetic allometry relations in newt populations. Nevertheless, it became apparent that in European 
newts, as well as in American newts (genus Notophtalmus, Harris 1989) small intrapopulation differences 
between transforming and non-transforming individuals in ontogenetic paths of size and shape chang- 
es developed, if any, and that happened relatively late in ontogeny. 


Acknowledgements 


We appreciate the help of Georg Dzukic, Davor Bejakovic and Ivan Aleksic in collecting specimens, as well as the 
help of Nikola Tucic in statistical analyses. 


References 


Alberch P., Gould, S. J., Oster, G. F. & D. B. Wake 1979. Size and shape in ontogeny and phylogeny. - Paleobiology 
5: 296-317 

Bookstein, F., Chernoff, B., Elder, R., Humphries, J., Smith G. & R. Stauss 1985. Morphometrics in Evolutionary 
Biology. - Academy of Natural Sciences, Philadelphia 

Cock, A. G. 1966. Genetical aspects of metrical growth and form in animals. - Quart. Rev. Biol. 41: 131-190 

Dzukic, G., Kalezic, M. L., Tvrtkovic, N. & A. Djorovic 1990. An overview of the occurrence of paedomorphosis in 
Yugoslav newt (Triturus, Salamandridae) populations. - Brit. Herpet. Soc. Bull. 34: 16-22 

Francillon-Vieillot, H., Arntzen, J. W. & J. Geraudie 1990. Age, Growth and Longevity of Sympatric Triturus cristatus, 
T. marmoratus and their Hybrids (Amphibia, Urodela): A Skeletochronological Comparison. - J. Herpetol. 24: 
1322 

Gould, S. J. 1977. Ontogeny and Phylogeny. - Belknap Press, Cambridge, Mass. 

Harris, R. N. 1989. Ontogenetic Changes in Size and Shape of the Facultativly Paedomorphic Salamander Notopiitial- 
mus viridescens dorsalis. - Copeia 1989(1): 35-42 


324 


Jolicoeur, P. 1963. The Multivariate Generalization of the Allometry Equation. - Biometrics 19: 497-499 

Kalezic, M. L., Crnobrnja, J., Djorovic, A. & G. Dzukic 1992. Sexual size difference in Triturus newts: geographical 
variation in Yugoslav populations. ~ Alytes 10: 63-80 

—- , Cvetkovic, D., Djorovic, A. & G. Dzukic 1994. Paedomorphosis and diffreneces in life-history traits of two 
neighbouring crested newt (Triturus carnifex) populations. - Herp. Journ. 4: 151-158 

Lande, R. 1985. Genetic and evolutionary aspects of allometry. In: W. L. Junger (ed.): Size and Scaling in Primate 
Biology, pp. 21-32. - Plenum Press, New York 

Lee, J. C. 1982. Accuracy and precision in anuran morphometrics: Artifacts of preservation. - Syst. Zool. 31: 266-281 

McKinney, L. M. (ed.) 1988. Heterochrony in Evolution: A Multidisciplinary Approach. - Plenum Press, New York 

Radovanovic, M. 1951. A new race of the Alpine newt from Yugoslavia. - Brit. J. Herpet. 1: 93-97 

Reilly, S. M. 1987. Ontogeny of the hyobranchial apparatus in the salamanders Ambystoma talpoideum (Ambystoma- 
tidae) and Notophtalmus viridescens (Salamandridae): The ecological morphology of two neotenic strategies. -J. 
Morph. 191: 205-214 

Shaffer, H. B. 1984. Evolution in a paedomorphic lineage. II. Allometry and form in the Mexican ambystomatid 
salamanders. - Evolution 38: 1207-1218 

Shea, B. T. 1984. An allometric perspective on the morphological and evolutionary relationships between pygmy (Pan 
paniscus) and common (Pan troglodytes) chimpanzees. In: R. L. Susman (ed.): The Pygmy Chimpanzee: Evolu- 
tionary Biology and behavior, pp. 89-130. - Plenum Press, New York 

—--— 1985. Bivariate and Multivariate Growth Allometry: Statistical and Biological Considerations. -J. Zool. Lond. 
(A), V 206: 367-390 

Tucic, N., Kalezic, M. L. & G. Dzukic 1985. Morphometric variability in a Triturus vulgaris population with Facultative 
Paedomorphosis (Amphibia). - Zool. Anz. 215: 102-108 j 

Verrell, P. A. 1985. Getting into a pickle with preserved specimens: formalin and distorsion in the smooth newt, 
Triturus vulgaris. - Herp. J. 1: 39-40 

Zar, H. J. 1984. Biostatistical Analysis. - Prentice-Hall, New Jersey 


Appendix 


Tab. 1. Mean values and standard error (X+SE) for analyzed forms of two smooth newt populations. n: sample size. 
See text for character and population designation. 


V1 V2 
metamorphic paedomorphic metamorphic 
larvae efts Ga 3d a8 őő larvae efts oho) od 
(n=53)  (n=48) (m=52) (n=52) (m=52) (n=52) (n=27) (n=34) (n=25) (n=25) 
X+SE X+SE X+SE X+SE X+SE X+SE X+SE X+SE X+SE X+SE 
iG, 36.85 35.63 80.56 79.42 75.92 72.40 PHAM 27.85 72.00 76.24 
+1.85 +1.88 +4.94 +4.34 +4.25 +3.88 +2.31 +2.02 +4.96 +5.48 
Lep 17.80 17.64 39.56 36.70 36.22 32.45 13.48 14.25 35.29 34.28 
+0.88 +0.73 +2.09 S715) +1.97 +1.56 +1.14 +0.84 asp l +2.38 
Led 19.05 17.98 41.00 42.72 39.71 39.95 13.63 13.60 36.71 41.92 
+1.34 +1.38 39.10 T299 +2.66 +2.85 +1.53 +1.35 E37 +3.69 
Lsv 20.00 19.78 43.29 41.18 39.88 36.96 15.01 15.91 38.96 39,15 
+1.09 +0.84 +2.41 +2.15 +2.08 +1.61 +1.22 +0.92 +2.66 +2.77 
Lte 3.99 3.74 7.60 7.04 7.18 6.62 3.08 3:02 5.91 5.92 
+0.18 +0.13 +0.45 +0.46 +0.45 +0.34 +0.25 +0.14 +0.39 +0.42 
Le 3.03 322. 5.56 577 4.81 4.89 2.53 2.74 4.96 5.20 
+0.21 +0.17 +0.31 +0.33 +0.30 +0.26 +0.24 +0.19 +0.31 +0.39 
Pa 6.22 6.48 14.46 16.03 1272 13.92 4.85 5.08 12.07 14.08 
+0.37 +0.35 +0.74 +1.26 +1.04 T095 +0.44 +0.43 +0.90 EISA 
Pp 5.72 6.27 14.05 16.10 1273 13.95 4.64 4.93 11.87 14.66 
+0.32 +0.33 +0.65 +1.24 +0.81 +1.01 +0.51 +0.33 +0.94 +1.41 
D 8.55 8.75 24.47 22.79 22.05 19.26 6.67 7.19 21.82 20.46 
+0.55 +0.82 +154 +1.30 +1.34 +1.28 +0.75 +0.49 L73 E7 
Lh 4.99 1.88 8.56 10.17 9.47 9.04 3.07 Sil 5.46 9.73 
+0.45 +0.11 IE) +1.47 +0.94 +1.00 +0.71 +0.14 +0.70 +1.55 


Tab. 2. Mean values and standard error (X+SE) for analyzed forms of two alpine newt populations. n-sample size. 
See text for character and population designation. 


Al A2 
metamorphic paedomorphic metamorphic 
efts ee 3d oe od larvae efts Ye) od 

(n=26) (n=25) (n=34) (m=54) (n=34) (n=50) (n=54) @=27)n=35) 

X+SE X+SE X+SE X+SE X+SE X+SE X+SE X+SE X+SE 

IL, 37.54 87.92 74.59 82.00 76.03 32.36 31.76 90.67 Wess 
+3.34 +4.90 +4.56 +4.60 +5.04 +1.85 +1.68 +5.03 +4.10 

Lep 19.96 45.53 38.73 42.74 38.79 16.95 17.08 45.83 38.17 
+1.11 +2.95 +2.45 +2.57 +2.74 +1.00 +0.88 +245 +1.97 

Led 17.58 42.39 35.86 39.26 3723 15.41 14.68 44.83 391.57 
+2.68 +3.05 +2.58 +2.84 +3.00 +1.30 +1.00 +3.37 +2.95 

Lsv 22.10 50.66 44.16 47.46 44.69 18.81 18.89 50.21 44.34 
+1.17 +2.81 +2.88 +2.98 +2.85 +1.04 +0.87 +2.42 +2.42 

Ete 4.49 9.10 8.19 10.13 9.86 4.07 3.89 8.57 7.74 
+0.19 +0.73 +0.65 +0.85 +1.16 +0.21 +0.19 +0.70 +0.52 

ite 3.47 6.36 6.05 6.21 6.03 2.89 3.11 6.13 5.89 
+0.22 +0.33 +0.44 +0.56 +0.51 +0.24 +0.21 +0.38 +0.29 

Pa Weave 16.78 16.54 FSi 1753 6.33 6.29 16.68 15.67 
+0.46 +1.89 +1.16 +1.69 +1.28 +0.43 +0.40 +0.82 +1.05 

Pp 6.98 17.03 16.31 16.64 17.14 5.75 5.85 16.31 15.41 
+0.44 +1.58 +1.28 +1.26 +1.25 +0.42 +0.46 +0.87 +0.89 

D 9.95 25.38 22.78 23.39 21.88 8.61 8.84 27-20 22.24 
+0.73 +2.91 +1.86 +1.88 +1.80 +0.65 +0.53 +2.26 +1.40 

Lh 2.15 6.29 6.32 AS 8.13 4.85 1.81 6.18 6.60 
+0.21 +1.53 +1.25 +1.74 +1.45 +0.44 +0.11 +1.22 +0.82 


Tab. 3. Mean values and standard error (X+SE) for analyzed forms of two crested newt populations. n: sample size. 
See text for character and population designation. 


(Gil C2 
larvae efts 2° 3d larvae efts 29 3d 
(n=35) (n=43) (=21) (n=19) (n=40) (n=90) (n=29) (n=20) 
X+SE X+SE X+SE X+SE X+SE X+SE X+SE X+SE 
E 83.89 89.67 130.81 123.89 64.18 73.93 142.48 122.80 
+3.53 +6.16 +9.31 +8.16 +5.50 +5.86 =+12.50 +10.40 
Lep 43.00 47.19 67.71 64.43 34.63 39.98 73.80 65.26 
+1.92 FAI +5.40 +4.36 +2.50 SRAM IAS) +6.16 +3.83 
Led 40.89 42.49 63.10 59.47 29.54 33.93 68.69 57.54 
E22 +4.28 +4.69 +4.41 +3.82 +3.92 +7,.00 +7.17 
Lsv 46.69 SESI 74.43 Z297 37.14 43.72 *81.66 73.84 
+2.05 +2.43 +6.20 +4.50 +2.74 +2.94 +7.01 +4.50 
Ete 9.19 8.77 12.84 12.46 7.39 7.34 14.00 12.49 
+0.53 +0.53 lel +0.96 +0.47 +0.56 aril oi +0.75 
Ike 5.91 6.26 8.30 8.09 5.05 5.54 8.79 TESI 
+0.31 +0.39 +0.68 +0.52 +0.27 +0.51 +0.73 +0.67 
Pa 16.41 17.36 24.43 25.79 12.82 14.85 252 25.70 
+0.80 +0.81 +1.88 +2.02 +0.94 +1.10 +1.78 +2.38 
Pp 16.76 17.95 24.98 26.49 12.69 14.35 25.50 26.39 
+0.77 +0.85 HS3 +1.67 +0.89 +1.12 +1.87 an) 38) 
D 23.44 26.89 40.56 38.71 18.23 23.04 44,93 39.97 
+1.68 +1.76 +4.31 +2.83 +1.65 +2.17 +4.52 +3.59 
Lh 8.35 722: 9.57 9.68 6.74 5.70 11.29 10.18 
+0.83 +0.62 +1.14 +1.36 +0.79 +0.70 +1.80 +1.14 


